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Abstract 

We consider a one-dimensional partial differential equation system mod- 
eling heat flow around a ring. The system includes a Klein-Gordon wave 
equation for a field satisfying spatial periodic boundary conditions, as well as 
Ornstein-Uhlenbeck stochastic differential equations with finite rank dissipation 
and stochastic driving terms modeling heat baths. 

There is an energy flow around the ring. In the case of a linear field with 
different (fixed) bath temperatures, the energy flow can persist even when the 
interaction with the baths is turned off. A simple example is given. 

PACS numbers: 44.10.-fi, 05.70.Ln, 05.10.Gg 
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Key words: Non-equilibrium statistical mechanics, stationary states. 

1 Introduction 

We consider the following system of coupled partial differential equations: 

dt(f){x,t) = T!-{x,t) 

dt7i{x,t) = {dl — l)(l){x,t) — g((f)(x,t)) — r]r(t) ■ a{x) 

dri{t) = -iri{t)~ri{ai,7r{t)))dt+y^idcOiit)t = 1,2. (1) 

In these equations is a field and vr is the corresponding momentum field; both 
satisfy periodic boundary conditions with x G [0,27r]. The non-linear term g is 
assumed to be bounded Lipschitz. The functions (actually distributions) a = (ai, 0:2) 
are fixed, real, and assumed to have Fourier coefficients di{n) oc for some — 1/2 < 
6 < 1/4, made precise below. (The brackets (■, ■) stand for the L'^[0, 27r]-inner product 
or distributional-test function pairing.) The vector r = (ri,r2) G is an artifact of 
the baths, and rj is a. coupling constant controlling the strength of the coupling of the 
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baths to the field. Finally, u = {uji,uj2) is 2-dimensional standard Brownian motion, 
and Ti and T2 are the bath temperatures. The system of equations ([T]) serves as a 
model of heat conduction around the ring. 

We will show, for the linear setting with g = and with the system in a stationary 
non-equilibrium state, i.e. with the temperatures Ti and T2 differing, that there 
can be a stationary energy fiow-a current-around the ring and, moreover, that this 
current persists even as the bath couplings are turned off, — ?■ 0. This phenomenon 
was noted but not proved in a previous article [18] with Y. Wang. The goals of the 
present article are to provide a proof of the claim and to give an explicit formula for 
the current. Parts of the analysis can be carried out assuming a non-linearity g in 
the wave equation, and so the non-linearity is retained in the analysis to the extent 
possible, with a view toward understanding its effect on this persistent current. 

The model itself has its origins in the heat flow models studied by Eckmann, Pillet, 
and Rey-Bellet [51 E] , who considered a flnite chain of non-linear oscillators coupled 
at either end to free fields that model heat baths. Assuming a judicious choice of 
couplings and Gaussian-distributed random initial conditions on the baths, these 
authors showed that, after integrating out the bath field equations, one is left with 
Markovian stochastic equations. They established existence of a unique stationary 
state exhibiting heat fiow and entropy production for these equations. See also [21 
El in [in 1151 US] for further developments of the model. The equations above, ([T]), 
are an adaptation of the models to an infinite-dimensional setting where the chain of 
oscillators has been replaced by a wave equation. 

Even in the case where g is cubic, g{y) = and thus unbounded, and the a's 
are chosen appropriately, the equations of motion ([T|) are known to have global in 
time solutions in spaces of low regularity, that is, with in H"', for any 7 > 1/3 

(7 = 1 corresponding to the usual energy norm) [17] . For the linear problem g = 0, 
the equations are of Ornstein-Uhlenbeck form. They have a unique invariant Gaussian 
measure (regardless of whether the temperatures are equal or not), and this measure 
is supported on fields (0, vr) such that (f){x) is almost surely Holder continuous with 
index < 1/2 [TS]. For bounded, Lipschitz g, the equations of motion with ultra-violet 
cutoffs have a unique invariant measure for each cutoff. By a soft argument, these 
measures have a weak* limit as the cutoff is removed, but it is not known whether 
the limiting measure(s) is invariant with respect to the above equations except in 
the equilibrium case Ti = T2, where a Gibbs state for the field is readily seen to be 
invariant [HI EOl [HI HH] • 

We turn now to the principal result. Let ai{n) be the n^^ Fourier coefficient of ctj, 

ai{n) = / e"*"^aj(x) dx. (2) 
Jo 

Assumptions. We assume that there exist positive constants Ci and C2 and another 
constant 6, — 1/2 < < 1/4 such that for all n 7^ 0, 

ci\nf <\ai{n)\<C2\nf i = l,2 (3) 
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and 

ci(|«i(n)|2 + |«2Hn < \al{n) + al{n)\\ (4) 

By way of explanation, the equations of motion are such that high frequency 
Fourier modes of the field couple rather weakly to the noise; however, the dissipation of 
these modes also is weak. Eq.Q insures that the modes couple sufficiently strongly to 
the dissipation to counterbalance the driving terms, thereby maintaining stationarity. 
(The situation is to be contrasted with stochastic heat equations and equations of 
hydrodynamics with viscosity, where high frequency modes dissipate very rapidly 
[U ttU 13 El IS]-) The second condition, Eq.(|l]), insures that degeneracy of mode 
frequencies is broken with the coupling a's turned on. As will be seen, the difference 
of frequencies between modes is the primary factor in controlling the correlations 
between field modes. 

The classical energy fiow through x G [0,27r] is given by TT{x)dx(f){x). We will 
consider this energy fiow averaged over the ring, i.e., J TT{x)dx(j){x) dx. A non- 
zero expectation of this quantity in a stationary state signals a net circular current. 
We project the equations of motion Eq.([T]) onto the subspace spanned by {e™^}, 
\n\ < M, and we let Em denote expectation with respect to the invariant measure 
for the system in this subspace. Although it is in fact a long calculation, we state the 
principal result ClS db theorem: 

Theorem 1.1 For the linear system g = 0, we have 
r 1 f'^'' 

lim lim Em — / 'n'{x)dx(p{x) dx 
7?-5>OM-!>oo 2n Jq 



(5) 



(Ti - T2) 1 >^ {al\n)al{n)) 

2 vr + l)|«?(n) + ajin^ + (|AfH| + lajinWY 

The theorem requires explanation. First, we emphasize that the result concerns 
only the linear problem g = 0. Second, for the equilibrium state or for a non- 
equilibrium limiting stationary state, 7r{x) is in H''~^, in H'', with 7 < 1/2 
almost surely, so that the pointwise heat fiow TT{x)dx4>{x) makes no sense; however, 
the product does make sense for the equations with an ultraviolet cutoff M. There 
is persistent current fiow for the finite-dimensional ultraviolet cutoff systems. The 
principal result here is the uniform convergence of this current, M — )■ 00 in 77 — )■ 0, 
and indeed that the current has a limit. It is this uniformity that we investigate 
here in detail, while including the effect of non-linearity where possible. Finally, we 
remark that when the coupling functions are turned off, the linear equations of motion 
certainly are not ergodic, and indeed invariant measures other than GIbbs states are 
possible. Of particular interest is how different non-equilibrium states are singled out 
in this limit, 77 — )■ 0. 
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The existence of such a current was announced in joint work with Y. Wang [T8] . 
What is new here is an account of the analysis behind the computation of the current, 
as well as a more explicit form for it in terms of the Fourier series coefficients of the 
coupling functions. We also provide an example below in which the coupling a's are 
5-functions at two points xq, xi G [0, 27r)- pictorially the situation in which the heat 
reservoirs are attached to a vibrating ring at just two points. The persistent current 
in this case shows a surprising discontinuous behavior as a function of the separation 
Xi — Xq of the two contact points. 

The result here, for a linear system, is in the spirit of the 1966 paper by Rieder, 
Lieb, and Lebowitz [TU] on stationary states for harmonic crystals. These authors 
considered a stationary state for a linear chain of oscillators with thermal reservoirs 
at either end maintained in general at different temperatures. They found that the 
expected kinetic energy of the j*'' oscillator is essentially constant in j in the interior of 
the chain (at the mean temperature), but that it dips as a function of j near the hotter 
end, with then a large discontinuous positive spike at the end oscillator. Similiarly, 
the expected kinetic energy grows in j in a neighborhood of the cooler end with a large 
discontinuous dip at the end oscillator. They also determined the expected energy 
flow, which is constant throughout the chain but which exhibits anomalous behavior 
when the coupling is large. Presumably these effects disappear with non-linearity, as 
might well happen to the current in our model if the non-linearity effectively impedes 
energy flow around the ring. 

We remind the reader that the existence of a stationary state invariant under 
the evolution equations for the non-linear {g ^ 0), nonequilibrium system remains 
an open question. One might attempt to construct an ultraviolet limit of stationary 
measures for g 0, M oo, in such a way that the expected current is fixed and 
finite. This would suggest how g is to be renormalized (if indeed renormalization is 
necessary) in order to achieve this limit. 



1.1 Example 

We conclude this introduction with a simple example in which the coupling distribu- 
tions, the a's, are chosen to be 5-functions, in particular 

ai{x)=6{x), 02(0;) = 0(^2 (x — xi) (6) 

with xi chosen < xi < vr, and c chosen < c < 1. In this case ai{n) = 1 and 
&2{n) = ce~*"^% and the growth (with ^ = 0) and non- degeneracy conditions on the 
a's, Eqs.( l3|^ . are seen to be satisfied. The limiting current is then given by 

r( ) = ^1 ~ ^2 'sp nc^ sin(2nxi) 

^^^>- 2tt ^ (n2 + !)(! + c4 + 2c2cos(2nxi)) + (l + c2)2' 

n:n>l ^ ' ^ \ / 

which is a conditionally convergent series. When the (5-functions are localized at 
diametrically opposite points, i.e., at and at xi = vr, this expected limiting current 
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C{xi) is zero, as one expects just for symmetry reasons. It is perhaps strange that 
C{xi) also vanishes for xi = tt/2. As seen numerically, C{xi) exhibits discontinuous 
jumps at both 7r/2 and vr but is otherwise continuous in xi. We have no physical 
explanation for this behavior. 



2 Calculation of the expected current 



2.1 Notation and eigenfunctions 

Let 



/ 



\ 



1 
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In terms of A(?7), the linear equations of motion Eq.([T]) read 

d^{t) = A(7/)$(t) + Vfdoo, 



(9) 



with the four-vector = t), 7r(a;, t), ri(t), r2(t))*. 

We will utilize the eigenfunctions and eigenvalues of A{r]) obtained via perturba- 
tion theory, A{ri) being a rank 2 perturbation of A(0). (See particularly the appendix 
of [19] for an account of the estimation of these eigenfunctions and their eigenvalues.) 
Note that A(?7) is not normal for rj ^ 0. We denote the right (column) eigenvectors of 
A(?7) by en{ri) = (en(/,(?7), en7r(^7), (??))*, the last component being an abbreviation 
for the two r components. The subscript n is a triple, n = (±,n, a), n > —1 being 
an integer, a = 1,2. The eigenvalue Xn{v) corresponding to en{r]) is equal to 



K{v) = ^±,nAv) = ±^(n2 + i)'/' 



2(±z(n2 + 1)1/2 + 1) 



(10) 



for n large (where by abuse of notation, 0{rin^) = 0{T])0{n^), for example). The 
eigenvalue shifts fin,a, cr = 1,2 are the eigenvalues of the 2 x 2-matrix 



M(n) 



\al{n)\ + \al{n)\ 
af{n) + d2^(n) 



df (n) + a2{n) 
\al{n)\ + \al{n) 



(11) 



and are given by 



(12) 
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For n = —1, the two corresponding eigenvalues of A(?7) are simply X-i^a = ~^+C>{r]'^); 
the ± plays no role. For the unperturbed rj = eigenfunctions, we take, for n > 0, 

e±,n,i(0) = -^(^^(n^ + l)"^/^cos(nx + ^„),cos(ra + ^),0)* 

e±,„,2(0) = -^(^^(n^ + l)-^/2sin(nx + Vn),sin(nx + V),0)*, (13) 
2a/ vr 



where ipn is the phase defined by 



al{n) + (n) 
\al{n) + d|(n) 



(14) 



For n = 0, e±fi = 2^(=F^5 1, 0)*; these eigenfunctions are not degenerate and there is 
no a index. For n = —1, e_i,i = (0, 0, 1, 0) and e_i,2 = (0, 0, 0, 1), the ± index is not 
in play. The momentum components of these functions, projected into the subspace 
of L^(0,2vr) spanned by e^*"^, are eigenvectors of M(n) above. 

Similarly, we have left eigenvectors (linear functionals) {/n('7)} of A{ri), their 
unperturbed versions taking the form 

/±,„,i(0) = -^{±i{n^ + lY/^cos{nx + ^n),cos{nx + 4j),0y 

j^,nM = -^(±2(n' + 1)^/2 sin(nx + V^„),sin(nx + V^),0)*, (15) 

with analogous expressions for /-t o and /-i,cr- Both the en(?7)'s and the /n(^)'s are 
normalized so that the real inner product {fn{v)j^n{i])) = 1 and so that, for n > 0, 
^n,n{v) has L^(0, 2vr)-norm equal to 1/2. We emphasize that we are particularly 
concerned with the large n behavior of these functions. Also, here and throughout, 
we are regarding these f^s as linear functionals so that, e.g., 

^iUiv)) = / {UAv,^)4>{x) + /n,,r(^?,a;)vr(x))c/x + UAv) " r. (16) 

We have that the momentum components of en(?7) and e°(0) differ by a small 
amount: 

Lemma 2.1 (Cf. [19j) For n = (±, n, a) 

en,n{v) - en,7r(0) = ^ an,k(^?)ek,7r(0) (17) 
k 

where the Fourier coefficients an,k('7) ore estimated 

0{T]^n'^^^'^ In^ n) k = n or k = (=l=, n, a) 
anM = { OivV'Hnn) k = (±, n, a'), or k = (t, n, a'), ^ (18) 
Oiv'M^) k ={±,k,a'), with k^n. 
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Moreover, we also have that the imaginary parts of these coefficients satisfy 

{Oirj^n^^-^ln^ n) k = n or k = (=F, n, a) 
0{rjV-Hnn) k= {±,n,a') or k= {T,n,a'), a' a (19) 
Ol'^^) k= i±,k,a'), with k^n. 

These estimates are uniform with respect to an ultraviolet cutoff M. 

We refer particularly to the appendix of [19j for details of these computations; 
although tedious, they just involve finite rank degenerate perturbation theory. The 
basic idea is to write the eigenvalue equations as two-dimensional implicit eigenvalue 
equations. We add: 

Remark. In these estimates and throughout this article one encounters sums of the 
sort ^fc>o A;^n fF^^- ^ Convenient method for estimating these sums is to break the 
sum up into three pieces: 1) an "infrared" part with k < n/2 and where the denomi- 
nator is simply replaced by the /c-sum then giving 2) a "near" part 
where n/2 < k < 3n/2, where the numerator is treated as a constant 3n^/2 and the de- 
nominator is replaced by n\n — k\ with the resulting J2k-n/2<k<3n/2 ~ 0(lnn), so 
that the near part then gives 0{n^~^ Inn); and 3) an "ultraviolet" part with k > 3n/2, 
where the summand is estimated simply by n^~^ with resulting sum 0{n^~^). Thus 
this example would give, overall, 0{n^~^\nn), furnishing some explanation for the 
mysterious logarithms in our estimates. 



2.2 Averaged expected energy flow 

In this section, we give an expression for the expectation with respect to the invariant 
measure of the energy flow 7T{x)dx4>{x) through a point x G [0, 2tt] averaged over the 
circle. We assume an ultraviolet cut-off M; Fourier modes of the field with frequency 
\n\ > M are simply set to zero, but all of our estimates are uniform with respect to 
the cutoffs, so the cutoff M is not explicitly expressed. Thus we consider 



E 



2tt 



n{x)dx(f){x)dx 



m n 



(20) 

where on the right side we have employed eigenfunction expansions for tt{x) and (pix) 
using the eigenf unctions of A(?7). The inner product is simply the L^(0, 27r)-inner 
product involving the vr and (f) components of the eigenfunctions. In order to prove 
the theorem, we will show that the diagonal terms, i.e. terms with m = n in this 
double sum, do not contribute in the limit 77 — )■ 0, nor do terms m, n with m ^ n. 
Only the near-resonant terms, terms such as m = (±,?t,, 1) and n = (±,n, 2), con- 
tribute in the limit, taking the form stated in the theorem. 
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Diagonal terms 



We begin by considering the diagonal terms in the above equation ( l20l) . terms with 
n = m. The estimates Eqs. fll8|19l) above and the fact that €^,^{1]) = Y^^n,n{v) from 
its eigenvalue equation imply that 

2i 

(en,7r(?7),C?xen,</,(r?)) = T-(3f?en,,r(^?),<9x5>en,,r(^?)) 

= Oir]W-Hnn) + Oir]^n^'-^). (21) 

Now in the linear theory, £'[|$(/n)P] is bounded by the supremum of the temper- 
atures, uniformly in r] small, see [IHl CH], so that {en,TT{v)y9xen,<t>{v))E[\^{fn)\'^] is 
certainly summable over n provided the exponent 6' < 1/4 as is assumed; the sum 
clearly vanishes in the limit — )■ 0. (In a non-linear theory, g 0, we would also need 
-£'[|'^'(/n)P] to be bounded or have suitably slow growth in n, but this is not known 
except in equilibrium, Ti = T2.) 



Off-diagonal terms 



We next consider the non-resonant off-diagonal terms m 7^ n with m ^ n ot the 
±'s differing in the double sum Eq.(l20l). The above estimates of the lemma, Eqs.fll7l 
[TSil . imply that 

{emAv),dxer.Av)) = oU^p^)+oU^^^^p^) (22) 

/ 1.29+1 \ 

The last term on the right side of this equation can be estimated using the Schwarz 
inequality (regarding one of the functions as fc^/(fc^ — m^)): 

''(^'™°""-') ^ E^^ O ( (F^S^) - (23) 

Thus 



(e„.,.(r/),9,e„,^(r7)) = O ^ 2 ) + ^ 



+0{7]^m''-W-'). (24) 

We also need an estimate on the correlation E[^{fm{v))*^{fn{v)]- Simply differ- 
entiating this expression with respect to time and using stationarity, we obtain the 
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identity 



'i^[<f(/m(r]))*(/„(r/),^7(0))] + i?[(/„.(r]),^7(0))*$(/„(r/))]) 

fL,M-TfnAv) 



(25) 



the last term in the third hne being an Ito-term. For this term fm,r{fl) ^-iid fn,r{v) 
are two-dimensional vectors, and T is to be interpreted clS cL 2 X 2-diagonal matrix, 
T = Diag(Ti,T2). For the linear theory, the first term on the right side is absent. 
The component fm,r{j]) = ''?(/m,7r(^?), a)/(Am(^?) + 1) is of order 0{rjn^~^); a similar 
estimate holds for fin,r{v)- We therefore have that 



for m 7^ n, at least in the situation where m ^ n. This quantity is certainly summable 
over m and n, m 7^ n using Young's inequality, the sum going to zero, 1] 0. 

The "antiresonant" terms in the double sum with m = n, but with different choices 
of the ± and no restriction on the a's, have matrix elements {em,TT{v) y 9xen,(f>{v)) 
which can be 0{1) in n. But the covariance factors are 0{rfin?'^~^), from the identity 
Eq.f E^ . with the denominator (A^(r7) + An(?7)) being 0{n). So again the products 
of the matrix element and covariance factors for these terms are summable, the sum 
going to zero, 77 — )■ 0. 

Remark: A term of the sort E[(^{f^{rj))* {fnij]), g{(j)))] appearing in the identity 
Eq.f l25p is no worse than K?0{n^^^'^)^ where k? is a uniform bound on variances 
i?[|$(/n(?7))p] and g is Lipschitz, a remark potentially of utility in analyzing the non- 
linear problem. 

These off-diagonal terms in the double sum also include terms with matrix ele- 
ments and field correlations involving the eigenvectors e_i^o-('7) and f-i^aiv)^ corre- 
sponding to eigenvalues near —1. In terms of the r-components of e_i^o-(^) that are 
(9(1) (actually nearly normalized), the 0-component is given by 




(26) 




(28) 
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The coefficient estimates Eq. (fT8|) then imply 

(en,7r(^?),<9xe-i,^,</,(?7)) = 0{rin^'^) 

(e-i,a,7r(^?),<9xen,0(?7)) = 0{r]n^''^) 

{e^i,^,M,d^e.i,^'Av)) = Oiif) (29) 
and, from the identity Eg. (125!) . 

EmuAvWHf-iym = o{i). (30) 

The resulting products such as {en,n{v)^ <9xe_i,o-,(/,(?7))-E [$(/n(77))*$(/_i,o-(?7))] are clearly 
summable in n and go to zero, r/ — )■ 0, hence do not contribute in the limit to the 
expected current. 

Near-resonant terms 

Finally we investigate the near-resonant terms in the double sum of the sort 

(e±,„,i,.(r/),a,e±,„,2,0(r/))S[$(/±,„,i(77))*$(/±,„,2(^))]. (31) 
We have by the Fourier coefficient estimates. Eg. (1181) . that 

{e±,n,iAv),dxe±,n,2Av)) 

(e±,n,l,7r(0),e±.„,l,,r(0)) + ( iv) (0))) 



•^±,ri,2('7) 

+ (e± ,n,l,7r iv) 

(0),e±,„,i,^(0)) ) 



,n,1.7r 

'^±,n,2{V) 



^ + Oiv'n''~Hn'n). (32) 

We have used 5a;e±^„_2,7r(0) = ne±^n,i,TT{0) and dxe±^n,i,n{0) = —ne±^n,2,TT{0)- Recall as 
well that the en(0)'s are normalized so that (e±,n,i,7r(0), e±,„,i_7r(0)) = 1/2, and that 
Xniv) = An(0) + 0(r/V^-i) with An(0) = ±i{n'^ '+ ly/^. We'summarize: 

T) 

{e±,n,iAv), dxe^,n,2Av)) = ± 2tin^ + iy/^ + Oiv'n''-') + Oiv'n''-' \n' n). (33) 
It remains to analyze the covariance factor £'[$(/-|- „ i(r7)*$(/-|- „ 2(^))] in Eg. fl3T]) . 



Again, we assume that there is no non-linearity in the identity Eg. fl25|) so that 

r±,nxM-Tf±,n,2Av) 



i^^[$(/±,n,i(r/))*$(/±,n,2(r]))] - 

(34) 



K,nAv) + ^±,nAv) 
V^{f±,n,lAv)^<^y ■ T{f±,n,2Av),<^) 



{K,nAv) + >^±,nAv)){K,nAv) + m±,nA{v) + 1) 



11 



{{fii,n{v)y(^i) is regarded as a two-component vector, i = 1,2). 

The numerator factor on the right side of this expression can be expanded: 

{f±,n,lAv),Oi)* ■T{f±,n,2Av),Oi) = -^^^^^ (/±,n,l,7r(^) , «)* " ^ (/±,n,2,7r(^?) , «) 

+ -^^^-^^(/±,n,l,^(^),«)* ■ (/±,n,2,^(^),«), (35) 

with D the two-dimensional diagonal matrix Diag(l, —1). But again, using the Fourier 
coefficient estimates Eq. fllSI) . which are also satisfied by the momentum components 
{/n,7r(^)} of the /n's, wc find this equal to 

^^^^(/±,n,i,.(0), a)* ■ /^(/±,„,2,.(0), a) + 0{iiV-' Inn). (36) 

Here, (/±,n,i,7r(^)5 and {f±,n,2,n{v)^(^) have nearly zero dot product (exactly zero 
for ?7 = 0), so that the (Ti + T2) part is represented in the small remainder terms. 
Moreover, via two-dimensional matrix algebra, 

(n\ ^\* n/f (n\ ^\ ^ (^f (^)«i(^)) /oy^ 

\ai{n) + ai{n)\ 

(The notation f„ is temporarily introduced to facilitate handling of already compli- 
cated expressions.) Thus the numerator factor is given by 

{f± 

(77), a) 

= ^^'~^'K n + 0{^V'Hnn). (38) 

The denominator of Eq. (1341) is expanded with the aid of the eigenvalue expansion 
Eq.([TO]): 



(A(±,n,i)(r/)* + A±,„,2(r/))(A±,„,,i(r7)* + l){X±^nM + 1) 
= 7 ^ , 7!.iu/2r , + 0{r,^n-Hnn). (39) 

Multiplying the factors ([38]) and ([39]) in ([MD and noting that z/„ = 0(^2^), we 
obtain 

(Ti - T2) 2zA^ 

2 ((/i„_i + /i„,2) ± ^(^^ + l)^/2(/i„,l - Ain,2)) 

+C(r/2n~2lnn). (40) 
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Equation ( 133|) and the previous Eq. (HO!) show that 

(e±,n,i,-(^),^xe±,„,2,<^(r7))^[$(/±,„,,i(r7))*$(/±,„,2(r/))] 
T1-T2 n Vn 



2 + 1)1/2 ((^2 + l)l/2(^^^^ _ ^^ 2) T ^Kl + /in,2)) 

+ 0(r/2n-2 Inn) + ^(r/V^-^) + ^(r/V-Mn^n). (41) 

In the hmit, r/ — )■ the following symmetries obtain: 

(e±,n,i,-(0),5x.e±,„,2,<^(0))E[$(/±,„,i(0))*<l>(/±,„,2(0))] 

(e±,„,2,.(0),9,e±,„,i,^(0))E[<l>(/±,„,2(0))*<f(/±,„,i(0))])* 

(e^,n,2,.(0), 9,e^,„,i,^(0))i?[$(/:p,„,2(0))*<f (/^,„,i(0))]) , (42) 



the first relation involving an integration by parts in the matrix element and using 
that en,</>(0) = j^en,7r(0). The third line follows from the facts that e±,„,a-,(/)(0) = 
—e^,n,a,,j>iS^) and an integration by parts so that 

(0), d^e^^nxM) = (e±,n,i,.(0), 9,e±,„,2,<^(0)), (43) 

and 

^[$(/T.n,2(0))*$(/^,„,l(0))] = E[$(/±,„,i(0))*<l>(/±,„,2(0))], (44) 

from Eq. fl34p . This simplifies the residual part of the double sum in Eq. fl2Up : the 
expected current is, in any case, real. 

In summary, the diagonal and off-diagonal terms and the remainder part of the 
near resonant terms in (HT]) in the double sum Eq. (l20|) are uniformly summable, and 
these sums go to zero, 77 — )■ 0. With the symmetry of the surviving terms, we have 
the final result 



lim E 

r)-s>0 



27r 

7r(x)(9j^0(x) dx 







27r 



2tx v^o 

m n 



(Ti - r2) 1 ^ Anl'n{^^n,l - /Wn,2) 



2 ((n2 + l)(/i„,i - Hn,2f + (/in.l + ^^n,2f) 



(Ti - T2) 1 >^ n'^{al\n)al{n)) 



2 vr ^f^^ (n^ + l)\al{n) + a^(n)P + (|«f(n)| + \al{n)\y 
by the definition of the /in,o-'s and z^n. This concludes the proof of the theorem. 
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